RNA-Seq-Based Transcriptome Analysis of Aflatoxigenic Aspergillus flavus in Response to Water Activity

Aspergillus flavus is one of the most important producers of carcinogenic aflatoxins in crops, and the effect of water activity (a(w)) on growth and aflatoxin production of A. flavus has been previously studied. Here we found the strains under 0.93 a(w) exhibited decreased conidiation and aflatoxin biosynthesis compared to that under 0.99 a(w). When RNA-Seq was used to delineate gene expression profile under different water activities, 23,320 non-redundant unigenes, with an average length of 1297 bp, were yielded. By database comparisons, 19,838 unigenes were matched well (e-value < 10⁻⁵) with known gene sequences, and another 6767 novel unigenes were obtained by comparison to the current genome annotation of A. flavus. Based on the RPKM equation, 5362 differentially expressed unigenes (with |log₂Ratio| ≥ 1) were identified between 0.99 a(w) and 0.93 a(w) treatments, including 3156 up-regulated and 2206 down-regulated unigenes, suggesting that A. flavus underwent an extensive transcriptome response during water activity variation. Furthermore, we found that the expression of 16 aflatoxin producing-related genes decreased obviously when water activity decreased, and the expression of 11 development-related genes increased after 0.99 a(w) treatment. Our data corroborate a model where water activity affects aflatoxin biosynthesis through increasing the expression of aflatoxin producing-related genes and regulating development-related genes.